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We present ab initio quasiparticle self-energy calculations in crystalline GaAs for cases of in- 
tense electronic excitation (~ 10% of valence electrons excited into conduction band), relevant 
for high- intensity ultra-short pulsed laser experiments. Calculations are performed using an out- 
of-equilibrium generalization of the GW approximation based on the Keldysh Green's function 
approach. Our results indicate that while the quasiparticle band gap is a sensitive function of the 
amount of excitation, it is not possible to induce complete band gap closure in this system by purely 
electronic means. 

PACS numbers: 78.55.Cr, 71.15.Qe, 71.10.-w, 71.55.Eq 

I. INTRODUCTION 

The ability to optically pump semiconductors with ultra-short pulsed lasers has led to the discovery of a host of 
phenomena over the past few decadesi. One of the most interesting is band gap renormalization (BGR), or a change in 
the quasiparticle gap of a material which occurs when electrons are redistributed in energy. The classic scenario is one 
in which electrons are excited from the top of the valence band to the bottom of the conduction band (in a direct-gap 
material such as GaAs), leaving holes behind. Though measurements of this effect are often difficult, reductions in 
the band gap have been inferred in the ~ 10 meV range for excited carrier densities of ~ 10^^ - 10^^ . Since 

these carrier densities are relatively small, the material remains crystalline during the duration of an optical pumping 
experiment. When the carriers have relaxed to the extrema of the bands, they are amenable to theoretical treatment 
within the effective mass approximation. The gap reduction itself is determined by calculating the quasiparticle 
self-energy in a screened- interaction (e.g. GW approximation^) perturbation theoretic treatment^'^. Extensions of 
these approaches which utilize ab initio electronic structure techniques have been applied to BGR in conditions of 
high doping levels, in which excess electrons (but no holes) are present. 

In the last few years, optical pump-probe experiments have been performed which reach carrier densities (electrons 
-I- holes) in excess of ~ 10^^ cm^^ ' these experiments, intense short-pulsed laser irradiation of GaAs rapidly melts 
the material and creates a metallic amorphous stated. Though the non-thermal melting is a primary goal of the 
studies, significant amounts of optical data were recorded prior to the time at which melting is thought to occur. For 
these short pump-probe delay times {'^ few hundred fs), deviations of the dielectric function from its values in the 
unpumped state are expected to be mostly electronic in origin (the heavier ions still being close to their equilibrium 
lattice positions). One of us (L.X.B.) has studied the changes in optical properties that result from the blocking of 
transitions (Pauli blocking) and the unbinding of exciton states which may occur in such conditions^. However in 
that work, BGR was neglected. 

Calculations of BGR in GaAs for these conditions of excitation have been performed by Kim et alt^A using an 
empirical pseudopotential technique together with a quasistatic effective mass-derived self-energy contribution. Their 
conclusion was that complete band gap closure would occur if 10% of the electrons are excited from the top of the 
valence band to the bottom of the conduction band. In their theory, a significant portion of the band gap reduction 
is due to the change in the screened electron-ion potential that results from the excitation of free carriers, the rest 
being caused by the changes in the quasiparticle self-energy. The change in Hartree potential (the interaction between 
electrons at the mean-field level) is neglected in their approach. 

In this work, we calculate BGR for laser-excited GaAs in which electrons have been excited from the top of 
the valence band to occupy the lowest energy states of the bottom of the conduction band. Single-particle states 
are taken from ah initio pseudopotential local density approximation (LDA) band structure calculations, and the 
quasiparticle self-energy is computed in the GW approximation (which includes dynamical screening of the electron- 
electron interaction)^. The LDA calculations are done for constrained (excited) occupation of the electronic states, so 
the self-consistent contribution of the Hartree term (owing to the redistribution of charge density) is included. We find, 
contrary to the results of Kim et al>iS, that it is impossible to induce complete band gap closure by purely electronic 
means, assuming that the excitation is of the type described above. The reason is that the change resulting from the 
electrostatic (electron-ion -I- Hartree) terms causes the gap to increase as electrons are moved from the valence band 
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top (near As) to the conduction band bottom (near Ga) . While our ab initio calculation of the self-energy contribution 
yields similar results to those of the more approximate theories2*&, it is the electrostatic contribution which dominates 
the picture at high carrier density (> 10%). In what follows, we describe the derivation of the self-energy contribution 
within the GW approximation using the Keldysh technique, present the details of our computation of BGR, and 
discuss our results. 



II. THEORY AND COMPUTATION 

A. Quasiparticle Self-energy Operator in the GW Approximation: Non-equilibrium Case 

We are interested in describing the single-particle excitations of a system of interacting electrons out of equilibriumii. 
Thus, we appeal to non-equilibrium quantum many-body theory, in particular, the Keldysh method. Using the Keldysh 
technique-^^.. Green's functions, self-energies, etc., for non-equilibrium problems can be constructed in a way similar 
to that for equilibrium problems. The only difference is the replacement of time-ordered Green's functions with 
contour-ordered Green's functions^S : 

G(ri, ra; n, ra) = -t{Tc[^Hiri;n)^+{^2;T2)]), (1) 

where the time-labels lie on a contour with two branches: t = t,ri {rj = 1, 2). (See Fig. ^) The contour C runs along 
the real axis, starting from -oo on branch 77 = 1, passing through ti and T2 once and returning to -go on branch rj = 2. 
Such a contour allows us to conduct the thought-experiment in which the system begins in some excited state of the 
non-interacting system, then the interactions are turned on adiabatically, turned slowly off again, and the system is 
returned to its initial state. This intellectual construct is necessary when relating expectation values of quantities in 
an interacting system to the corresponding expectation values of the non-interacting systemic. In what follows, we 
omit the spatial variables for simplicity, which are to be used exactly as in the equilibrium case. We return to the full 
notation later. 

The contour-ordered Green's function obeys the same Dyson equation as the time-ordered Green's function, if one 
replaces the real time axis integrals by contour integrals: 

G(Tl,T2) =Go(Tl,r2)+ / dT3 [ dT4Go(Tl,T3)S(T3,r4)G(T4,r2). (2) 

JC Jc 

We use the GW approximation for the self-energy: 

I](ti,T2) = zG(ri,T2)VK(ri,T2), (3) 
where the screened Coulomb potential W obeys the following self-consistent equation: 

W{n,T2) = V5{ti -t2) + V [ dT3P(ri,r3)l¥(r3,T2), (4) 

Jc 

where V is the Coulomb potential. The polarization function is calculated within the random phase approximation 
(RPA): 

P(Tl,T2) = -^G(Tl,T2)G(T2,Tl). (5) 

In order to make the calculations tractable, we replace the contour integrals by real time integrals. This is achieved 
using the Langreth rules for analytic continuation^^, which we present in the Appendix, after introducing some 
commonly used notation for functions. A, defined on the contour (called 'lesser', 'greater', 'retarded' and 'advanced' 
respectively) : 

A<it,t') = A{t,r^ = l,t\r^' = 2) 
A>{t,t') = A{t,rj^2,t\rj' = l) 
A''{t,t') = e{t-t')[A>{t,t') - A<{t,t')] 

A"(t,t') = e{t' - t)[A<{t,t') - A>{t,t')]. (6) 

Our goal is now to evaluate the retarded quasiparticle self-energy operator, E^, in the GW approximation. From (3), 
(6), and (A. 4) (see Appendix), S'' becomes: 

T.^'it, t') = i[G<{t, t')W''{t, t') + ^{t, t')W< {t, t') + G'{t, t')W'"{t, t')l 
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which can be further simpHfied using (A.l) and (A. 2): 

E'-(t,<') =i[G<{t,t')Wit,t') + G''{t,t')W>it,t')]. (7) 
The Keldysh components of the screened Coulomb interaction are found using the Langreth rule (A. 3): 



W'^it, t') = V5{t ~t') + V / dtiP'^it, ti)W{ti,t') 

J — oo 

f'OO 

W"(i, t') = V5{t -t') + V / dtiP^it, ti)W''{ti,t') 



W<{t,t')^V / dh[P^{t,h)W<{t[,t') + P<{t,h)W^{h,t')] 

J — OO 

W>{t,t')^V dh[P''{t,t,)W>{t[,t') + P>{t,h)W''{h,t')]. (8) 



It can be shown^^-^^ that the self-consistent equations (8) lead to the following relations for the "greater" and 
"lesser" components of the screened interaction: 



W<it,t')= dh dt2W{t,ti)P<iti,t2)W^{t2,t') 

J — OO J — OO 

poo poo 

W>{t,t')^ dh dhW'{t,h)P>{tut2)W''{h,t'). (9) 



We construct "non-interacting" Green's functions within the quasiparticle approximation from LDAiS. Bloch wave- 
bmctions |rik), and eigenvalues £'„k- Using (1) and (6) we find: 

(rik|G"'(tj)|mk'} = — -—SnmSkw 

- Aik + if] 

{nk|G"(w)|TOk') = J- -5,„„(5kk, 

(nk|G<(tj)|TOk') = 2TTif{Enk)S{0J - £^„k)'5nm<5kk' 

{nk\G>iLu)\mk') = -27rz[l - /(^„k)]<5(cc' - S„k)'5n™4k', (10) 

where i] = 0^ and we assumed that the "non-interacting" system can be characterized by a distribution function in 
energy, f{E). 

In frequency space, the expression for the retarded self-energy (7) becomes: 

(7'ik|S]'"(a;)|nk) = z AfQ(n, m, k, — q)MG' (n, m, k, — q) 

m,q,G,G' 

/CXD T 
^ (q, 0J)G< (k - q, £;„k ~oj) + W> ,,, (q, c.)G'^(k - q, E^y, - cu)], (11) 

where the matrix elements M are defined by: MG(".,TO,k, — q) = (mk — q|e~'^'i+*^-''"|nk). Next we need the Fourier 
components of the screened Coulomb interaction. Using (8) and (9), the Fourier components of the screened Coulomb 
interaction can be written: 

I¥4G'(q,^) = [e'-(q,c^)-']G,G'^^(q + G') 
W^G,G'(q,^) = [e"(q,c.)-']G.G'V^(q + G') 
W< G,(q, Lo) = -[e'-(q, Lo)-h< (q, c.)e'^(q, c.)-1]g.G' V^(q + G') 

W> ,,, (q, u;) = -[e'-(q, u;)-h> {q, c.)e'^(q, c.)-i]G,G'y(q + G'), (12) 
where the dielectric function is defined by: 

ec.G' (q, ^) = '^G,G' -vi(i + G)pa,G' (q, ^) 

eG,G' (q, ^) - -^G^G' - V^(q + G)P^,G' (q, ^) 

4,G'(q>^) = -nq + G)P<G.(q,^) 

e> ,,,(q,c^) = -y(q + G)P>G,(q,^). (13) 
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The Fourier components of the polarizabihty are found via (5), (A. 5) and (10): 



-PG,G'(q>^)= E ^^-G("."i,p,q)M_G'('^>™>p>q)- 



/(-B„p+q) - f{E.mp) 



Enp+q ~ LlI — Ernp + 
n,m,p r- I -1 t- 



-PG,G'(q'^) = -27ri E Ml(.{n,m,p,c[)M^G'{n,m,p,(i)f{Enp+q)[l-fiErap)]5{E,ap+q-uj-E„,p) 

n,m,p 

-PG,G'(q''^) = -27ri E A^-G(":"z,p,q)M_G'(n,m,p,q)[l-/(£;„p+q)]/(^,„p)5(£;„p+q-w-i;™p). (14) 



n,m,p 



Decomposing the self-energy into the screened-exchange and Coulomb-hole parta^: 

{n\i\i:''{E)\nk) = {nk\j:'sxiE)\nk) + (nk|I]j;^(£;)|nk), 
the final expression for the retarded self-energy is: 



(nk|I]5_Y(^)l'^k) = — E M^{n,m,'k, ~q)MG'{n,m,'k, ~q) 

m,q,G,G' 

x[e'^(q,^ - ^,„k-q)"']G,G''^^(q+ G')/(£;„.k-q) (15) 
{nk\T,cH{E)\nk} ^ - E ^^g("; k, -q)AfG'("i ™> k, -q) 



m,q,G,G 



/OO 7 -1 
-^^:^[e'-(q,c.)-le>(q,..)6'^(q,c.)-l]G.G'r;^ — ^q + G'). (16) 

The above result represents a general expression for the self-energy (within the GW approximation) of a system 
characterized by an arbitrary distribution function, f{E). In particular, in the case of a system in equilibrium at 
some temperature, f{E) is just the Fermi distribution function, fp, and the following holds via (13) and (14): 

eG.G'(q'^) = eG,G'(q.^)[i + tttttI' 

where /b denotes the Bose distribution function. Using this result together with the Keldysh relation (A.l) for e , it 
can be shown that (16) becomes: 

m,q,G,G' 

1 r ,Je-(q,..)-^-.-(q,^)-^]G,G> l + fsH _^ (^^^ 



X ■ 



TT 7_oo 2i S - - -Bmk-q + IT] 

For systems with inversion symmetry, we can further replace in (17) ^'^ — ~^ (q.'^) — 1 |-,y Jrfi[e^{q,uj)~^] and thus 
we recover a previous result obtained by applying the finite temperature Matsubara Green's function formalism to 
a system with inversion symmetry. 



B. Computational details: GaAs 

We calculate the quasiparticle (QP) band-structure for crystalline GaAs in various states of electronic excitation 
by first order perturbation theory^ from the real part of the self-energy: 

£;^Sr = + {nk\Rei:{E^^) - V^^lnk), (18) 

with Re's = ^(S'' -I- S°). Using {nk\ReT,\nk) — i?e(7ik|S''|nk), (18) is solved by evaluating the real part of the 
retarded self-energy, (15) and (16), for various distribution functions /. 

The "non- interacting" mean field wavefunctions |nk) and eigenvalues Enk as well as the exchange-correlation 
potential Vxc are obtained within constrained LDA performed at the experimental lattice constant. We use ab initio 
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pseudopotentials generated with the scheme of Trouiller and MartingSS. A plane wave basis with an energy cutoff of 
30 Ry is used to represent the Bloch states. We use 50 valence plus conduction bands and discrete Monkhorst-Pack 
meshesSi of q points (8x8x8) in the first BZ in the sums of (15) and (16), while crystalline local-field effects in the 
dielectric matrix were included by summing over G, G' up to a cutoff of 6 Ry. The integral over lo in (16) is evaluated 
by setting the integration limits to ±50 eV. 

We consider five occupation number distributions /, each corresponding to a different excited carrier density: 0%, 
5%, 10%, 15% and 20% of the total number of valence electrons in the system (the coarseness of our meshes of q points 
precludes us from considering excited carrier densities between 0% and 5%). Each is characterized by a quasi- Fermi 
level of the conduction bands: the energy below which all conduction states are filled. The energy above which all 
valence states are empty is then determined by requiring charge neutrality. Note that our assumptions imply that the 
individual excited electron and hole populations « are in quasiequilibrium at T=0. Fig. |21 shows the excited carrier 
occupation in the Brillouin zone (BZ), for the various excited carrier densities considered in this work. Because the 
effective mass of the conduction band at F is very small, the F valley of the conduction band can accommodate less 
than 1% of the total number of valence electrons (jitot)- It follows that when the excited carrier density is several 
percent of ntot, the excited electrons must reside in a combination of conduction F, L and X valleys. It is these high 
carrier densities that we are concerned with here, given the conditions reached in recent pump-probe experiments^. 

III. RESULTS AND DISCUSSION 

Before presenting the results of our calculations of the QP self-energy operator, we discuss the changes in the 
band gap caused by the excitation at the mean-field (in this case, LDA) level. Fig. 01 presents the constrained LDA 
calculation results of BGR. We see that the valence-conduction gap E^^^ increases linearly with the excited carrier 
density, the slope being only slightly smaller for the direct gap at F than for the indirect gaps at L and X. In order 
to understand these results, we analyze the dependence of the direct gap at F on the excited carrier density. We 
do this by decomposing the Kohn-Sham hamiltonian into the kinetic energy part K, the ionic potential plus the 
Hartree potential Vion + VHartree, and the LDA exchange-correlation potential V^c- Their individual contributions to 
BGR are presented in Fig. 0] which shows that the dominant term leading to the increase of the constrained LDA 
valence-conduction gap at F is the Hartree contribution. Why is the Hartree term in the band gap increasing with the 
excited carrier density? The answer is straightforward if we take into account that: 1) high- lying valence electrons 
are localized more near As, while low-lying conduction electrons are localized more near Ga. 2) under excitation the 
electron charge density goes from As to Ga. Then it follows that under excitation, the electrostatic energy between 
one electron and the total charge density goes down for valence electrons and up for conduction electrons, resulting 
in an increase in the valence-conduction band gap. We note that at the mean-field level, our results are different from 
those obtained by Kim et ali^ using a self-consistent screened empirical pseudopotential treatment. They found (for 
10% excited carrier density) a significant decrease of the indirect gaps 1 eV for the F-X gap) and no change of the 
direct gap at F. We ascribe this to their neglect of the change in the Hartree term. Though some of this mean-field 
effect is undoubtedly included through the screening of their pseudopotentials by the excited carriers, we suspect that 
a large portion of this Hartree interaction is incorrectly accounted for in their approach. 

We now move to results at the QP level, by including the self-energy corrections in the GW approximation. Fig. 
|S1 presents the main results of our work, showing that the BGR picture changes dramatically when we go from the 
mean-field level (constrained LDA) to the QP level: we have a substantial decrease of all band gaps for relatively low 
excited carrier density (less than 5%), then the band gaps are rather constant for excited carrier densities between 
5% and 10%, while for relatively high excited carrier densities (larger than 10%) the band gaps increase linearly. We 
note that the lowest value for the renormalized gap is about 0.8 eV, which results from roughly a 10% excited carrier 
density and is still a direct gap at F. 

We analyze the results of Fig. Elby separating out the contribution from the exchange-correlation self-energy part. 
We do that by introducing 5J:g{n%) = T,g{n%) - Sg(0%) (with Sg = {ck\ReT,\ck) - {vk\ReT.\vk) , v and c denoting 
the valence and conduction bands near the Fermi level, k being the F, L or X) and we interpret JSg as the self-energy 
contribution to BGR. Fig. shows the dependence of ST,g with the excited carrier density. We see that ST,g is 
negative and roughly independent of the gap we are looking at (it is largest for the indirect gap from F-X), and that 
it decreases rapidly with the excited carrier density for relatively low densities. For higher excited carrier densities, it 
becomes rather independent of density. It is worth noting that Kim et al. found a similar self-energy contribution for 
an excited carrier density of 10%: (5I]g(10%) « —0.9 eV (an average along the entire F-X direction), with the largest 
value at X. 

There are two main sources for the self-energy contribution to BGR. One comes from the change in the Fermion 
occupation numbers, /, appearing explicitly in eq. (15). The other comes from the change in the electronic screening 
which manifests in both eqs. (15) and (16) through the dielectric function e (of course, e depends implicitly on /). 
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By increasing the excited carrier density, the screening acquires a metalHc character due to the increased number of 
'free' electrons in the conduction band and 'free' holes in the valence band. We find that the change in the electronic 
screening (affecting both the screened-exchange and the Coulomb-hole terms) accounts only for a small fraction of 
the BGR. For the excited carrier densities considered in this work (in excess of 3.6 x 10^^ the maior source for 

the self-energy contribution to BGR was found to be the change in the Fermion occupation numbers which enter the 
screened-exchange self-energy expression (see the / in eq. [15]). This picture may not hold for excited carrier densities 
lower than 10^^ trn^ ^- illustrate this observation in Fig. [3 where we show the screened-exchange self-energy 
contribution to the direct gap at F calculated in two ways: the solid line represents self-consistent calculations based 
on eq. (15), while the dotted line represents calculations based also on eq. (15), but in which the dielectric function 

was computed for the unexcited case. Finally, we note that our results for ^S^"^ follow an approximate ^ power 
law behavior in the excited carrier density; the same power law behavior for ^I'^^g^ can be obtained employing a very 
simple model based on the parabolic band approximation and replacing e'' with a simple dielectric constan1iS*iS. 

Summarizing our results for BGR, we conclude that for relatively low excited carrier densities (up to 10% of the 
total valence electron density) the QP valence-conduction band gap decreases (but never closes) primarily due to the 
change in the Fermion occupation numbers entering the screened-exchange QP self-energy expression, while for higher 
excited densities the gap reaches a minimum value and undergoes an upturn with increasing density. The increase 
is mainly due to the change in the electrostatic energy at the mean-field level described by the Hartree term. Our 
prediction that there is no band gap closure due to these valence band top — > conduction band bottom excitations 
agrees with the measurements of L. Huang et alf2', in which the dielectric function of GaAs was measured at short 
100 fs) time delays of a pump-probe experiment and no signature of metallicity was found^'^. We also mention that 
the approximate k-independence of the BGR that we predict suggests that the overall shape of the dielectric function 
(i.e.- the distribution of oscillator strength in the optical absorption spectrum) in these excited configurations will be 
similar to that predicted by one of us^ using an approach in which BGR was neglected. Finally, we add that our 
computations of the imaginary part of the QP self-energy for the cases considered above yield values small enough 
for us to conclude that lifetime broadening due to the scattering of excited carriers should be of little importance in 
determining the shape of the optical spectra. 

IV. CONCLUSIONS 

We have studied band gap renormalization in laser-excited GaAs. The cases we considered were those of intense 
excitation (~ 10%) in which excited electrons occupy the lowest-lying conduction states and the excited holes occupy 
the highest-lying valence states. The quasiparticle self-energy was computed with a non-equilibrium variant of the 
GW approximation using the Keldysh technique. Our study indicates that it is not possible to induce complete band 
gap closure by purely electronic means. In reaching this conclusion, we found that the contribution of the Hartree 
term to BGR, describing the mean-field interaction between an excited electron (or hole) and the charge density of the 
remaining electrons, is important and cannot be neglected. Our findings seem to support experimental observations^ 
in which no evidence of electronically-induced band gap closure was seen. 
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Appendix 

In this section, we present the Langreth rulei^ for analytic continuation of products of functions integrated over 
contour C. This is necessary when evaluating the GW self-energy expression. The functions introduced in (6) satisfy 
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the Keldysh relation: 

A'^{t,t') ^ A''{t,t') = A>{t,t') - A<{t,t'), (A.l) 
and also that for any functions A and B: 

^'■(t, t') = A'^it, t')B''{t, t') = 0. (A.2) 

The Langreth rule for analytic continuation states that if on the contour: 

F ^ AB 
Jc 

then on the real time axis (for a concise derivation seei&): 

/oo 
A^'B'' 
-OO 

F" = / A'^B" 



F< = / A''B<+A<B'' 

' —OO 
nOO 



F> - / +yl>S^. (A.3) 



It can also be shown that if: 



then: 



and: 



F{t,t')^A{t,t')B{t,t') 
D{t,t')^A{t,t')B{t',t), 



F'-(i,t') = A<{t,t')B'-{t,t') + A'-{t,t')B<{t,t') + A'-{t,t')B'-{t,t'), (A.4) 



D''{t, t') = A< {t, t')B''{t', t) + ^'■(t, t')B< {t', t) 
D'^it, t') = A< {t, t')B''{t', t) + A'^it, t')B< {t', t) 
D<{t,t') = A<{t,t')B>{t',t) 

D>{t,t') = A>{t,t')B<{t\t). (A.5) 

This result can be directly applied to the evaluation of the GW self-energy expression, where G and W replace A and 
B respectively (as well as to the polarizability P, replacing A and B by G). 
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FIG. 1: Closed time path contour C along real axis, showing two branches, 77 = 1, 2. 
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FIG. 2: Excited carrier occupation in GaAs; the LDA band structure shown here is taken from the unexcited case. 




FIG. 3: Theoretical Kohn-Sham band gaps of GaAs calculated within the constrained LDA. 
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FIG. 5: The quasipaxticle band gaps of GaAs as a function of excited carrier density calculated within the GW approximation. 




FIG. 6: Self-energy contribution to the band gap of GaAs as a function of excited carrier density. 



